Super-sharp resonances in chaotic wave scattering 
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Wave scattering in chaotic systems can be characterized by its spectrum of resonances, z n = 
E n — where _E n is related to the energy and F n is the decay rate or width of the resonance. 
If the corresponding ray dynamics is chaotic, a gap is believed to develop in the large-energy limit: 
almost all T n become larger than some 7. However, rare cases with T < 7 may be present and 
actually dominate scattering events. We consider the statistical properties of these super-sharp 
resonances. We find that their number does not follow the fractal Weyl law conjectured for the 
bulk of the spectrum. We also test, for a simple model, the universal predictions of random matrix 
theory for density of states inside the gap and the hereby derived probability distribution of gap 
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Scattering of waves in complex media is a vast area 
of research, from oceanography and seismology through 
acoustics and optics, all the way to the probability am- 
plitude waves of quantum mechanics [lHj] . We focus our 
attention on the important class of systems for which 
the complexity is not due to the presence of random- 
ness or impurities, but rather because the correspond- 
ing ray dynamics is chaotic. The presence of multiple 
scattering leads to very complicated cross-sections, with 
strongly overlapping resonances that, although determin- 
istic, have apparently random positions and widths. It is 
not uncommon that in a given situation only the sharpest 
resonances are relevant, with the others providing an ap- 
proximately uniform background. 

For concreteness of terminology, we consider quantum 
mechanical systems, but our results are general. We de- 
note resonances by E — i~ and call T the width or the 
decay rate. In the large-energy limit a gap develops in the 
resonance spectrum: most resonances have their widths 
larger than 7, which is the average decay rate of the cor- 
responding classical (ray) dynamics. This was noticed 
long ago [5| , and more recently there have been attempts 
to prove it rigorously [gj]. We are interested in the rare 
case of states inside the gap, i.e. with r < 7, which we 
call super-sharp resonances. 

The distribution of typical resonances in chaotic sys- 
tems is conjectured to follow the so-called fractal Weyl 
law Q : their number grows with E as a power law, whose 
exponent is related to the fractal dimension d of the clas- 
sical repeller, the set of rays which remain trapped in 
the scattering region for infinite times, both in the fu- 
ture and in the past Q. In a numerical experiment with 
a simple model, we find that the number of super-sharp 
resonances, denote it by Nssr, does not follow this law. 
It does grow with E according to a power law, but the 
exponent seems to be insensitive to 7 or the dimension 
of the repeller. 

We also investigate the dependence with energy of the 
width of the sharpest (and usually most important) reso- 
nance. Let this be denoted Tq. As E grows, it is expected 



to converge to 7. More natural variables are their expo- 



nentials, and we find that e 



e 7 decreases with E 



according to a power law, whose exponent is well approx- 
imated by d — a. 

A very fruitful approach to chaotic scattering of waves 
is random matrix theory (RMT) in which the sys- 
tem's propagator (the Green's function of the wave equa- 
tion) is replaced by a random matrix, whose spectral 
properties are studied statistically. The RMT prediction 
for the density of resonance states inside the gap was 
derived in 10]. In this work we derive the RMT predic- 
tion for the probability distribution of e~ r ° — e~ 7 , and 
compare both these predictions to numerical results in a 
specific system. 

Waves in chaotic systems can be modeled by the so- 
called 'quantum maps'. These are N x N unitary matri- 
ces where N ~ and the large energy limit E — > 00 is 
replaced by the limit of large dimension, N — > 00. This 
approach has been used to study transport properties of 
semiconductor quantum dots I 111 , entanglement produc- 
tion Il2[, the fractal Weyl law [13j, fractal wave functions 
14, Il5j . proximity effects due to superconductors [l6j . 
among other phenomena. Scattering can be introduced 
by means of projectors. The propagator becomes a sub- 
unitary matrix of dimension M < N, and its spectrum 
comprises N — M zero eigenvalues and M resonances of 



the form z n 



i(E n -ir n /2) 



As our dynamical model we use the kicked rotator. Its 
classical dynamics is determined by canonical equations 
of motion in discrete time, 

K 

qt+i = qt+Pt+ ~r sin(27rg f ) (mod 1) (1) 

47T 

K 

Pt+i =Pt + t~ [sin(27rgt) + sin(27rg t+ i)] (mod 1). (2) 

47T 

This system is known to be fully chaotic for K > 7, 
with Lyapunov exponent A ~ ln(if/2). Its quantization 
H] yields a N dimensional unitary matrix U. We 
set equal to zero a fraction of 1 — \l of its columns, corre- 
sponding to a 'hole' in phase space. Here [i corresponds 
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FIG. 1: Scaling exponents as functions of /i for the open 
kicked rotator, d, a and ft are related to the fractal Weyl law, 
the number of super-sharp resonances and the width of the 
gap, respectively. We see that a is approximately indepen- 
dent of the size of the opening. The solid line is a best linear 
fit to d, while the dashed line is d — (a), where (a) ~ 0.66 is 
the average value of a. 



to the fraction of rays which escape the scattering re- 
gion per unit time, so fi — e~ 7 . We do not consider very 
small values of fx, which would correspond to widely open 
systems. 

We find numerically that AfssR. ~ N a , while e~ r ° — 
e ~T ~ . The exponents are plotted as functions of fx 
in Figure 1, together with d, the numerically determined 
exponent in the Weyl law. We see that, somewhat sur- 
prisingly, in this range a is approximately constant, i.e. 
insensitive to the dimension of the repeller. The super- 
sharp resonances do not follow the fractal Weyl law. On 
the other hand, the exponent j3 has approximately the 
same slope as d. We find the relation d — j3 = a to be 
approximately fulfilled, which is to be expected since it 
says that the number of super-sharp resonances, ~ N a , 
is proportional to the total number of resonances, ~ N d , 
times the width of the gap, ~ . 

We now turn to an RMT treatment of the problem. 
RMT for quantum maps amounts to taking matrices 
uniformly distributed in the unitary group. In (l7| . an 
ensemble of truncated unitary matrices was introduced 
as a model for scattering, and it was shown that, as 
N,M -> oo with n = M/N held fixed, the probability 
density of r = e~ r / 2 converges to 



Pco(r) = 



1-fx 



2r 



Li (1 - r 2 f 



(3) 



if r 2 < /J, and to otherwise. The expression ((3]) was 
tested numerically for a chaotic quantum map in (lij 
and found to be an accurate description of the bulk of 
the spectrum, provided a rescaling was performed to in- 
corporate the fractal Weyl law. 

Let us start with the density of super-sharp resonances. 



1.0 

0,8- 

0,6 

0,4 

0.2 

0,0 4. 






N=500 




N=1500 




N=5000 




N=10000 



0.0 



0,1 



0,2 



0,3 



0.4 



0.5 



FIG. 2: Top panel: Pi = Pr(r 2 — fi < x) for the open kicked 
rotator, at several dimensions, for fi = 0.8. Lower panel: same 
for P2 = Pr(r;(r 2 — fx) < x), where rf is defined in the text. 
Solid line is a rescaled RMT prediction. 



This calculation can be found in [lOj; we sketch it here 
for completeness. The probability density for r = e~ r / 2 
is known exactly [13] , and when M = fiN and N ^> 1 it 
can be approximated to 



P(r) 



Poo(r) 



1 + erf rj^/Jl 



.(fx-r 2 ) 



where erf is the error function and 



N 



2fi{l-fx) 



(4) 



(5) 



is now the large number. Clearly, (Q} will be different 
from Pqo (r) only if /1 — vq is of the order of 1 jrf. After 
setting r 2 = fi + e/rj, the probability distribution for the 
variable e becomes 



P(e) = v^F(l - erf(e)). 



(6) 



This is the density of states inside the gap. It is only 
non-zero in a small region that shrinks as N^ 1 / 2 in the 
asymptotic limit. Its integral provides the probability 
that r 2 — fx be less than some value x. This we denote 
by Pr(r 2 — /1 < x). 
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In Figure 2a we see Pr(r 2 — \x < x) for the open kicked 
rotator, at various dimensions for \x = 0.8. We see that 
as N grows the values of r become more localized around 
/j. In Figure 2b we introduce a scaled variable rj(r 2 — /i), 
but with rj different from r\ in that it involves the actual 
exponent /3 instead of the RMT prediction 1/2: 



V 



N 



2/x(l - a*) 



(7) 



All curves fall on top of each other, indicating that this 
is the correct scaling. Moreover, the shape of the curve 
agrees with ((6]). 

Let us now consider the probability distribution of the 
largest eigenvalue of the propagator, which corresponds 
to the sharpest resonance and whose modulus we de- 
note R. Largest eigenvalue distributions are important 
in different areas of mathemati cs Jl9| and physics [2fjj |. 
and have even been measured [21] . A similar calcula- 
tion to the one below can be found in (2^. Let V({z}) 
denote the joint probability density function for all eigen- 
values. When integrated over all variables from to x, it 
gives the probability that all eigenvalues are smaller than 
x. Therefore, if Q{x) is the probability that the largest 
eigenvalue be smaller than x, then 



Pr(R < x) 



M 



[0,x] 



V{{z})\{d 2 



(8) 



The jpdf of the eigenvalues is [Tj 



1..M 



n 



2\N-M-1 



(9) 



i<3 



It is a usual trick to write Tii<j \ z 



zj\ 2 in terms of the 

Vandermonde determinant |detA| 2 , where = z l f . 
This can be shown to be equal to MldetB, where Bij — 
zf~ (z*) 1-1 . Each element of the matrix B depends on 
a single variable, and the integration decouples. The an- 
gular part of the integrals diagonalizes the matrix and, if 
M = fiN and N ^> 1, the result becomes 



Px{R < x) cx TT / (1-y) 



y 3 dy. (10) 



This result is exact, but a bit cumbersome. Approximat- 
ing the integrand by a Gaussian function we arrive at 

Pr(i? < x) cx Utio {l + W [v - I* + ^)] } • 
This result can be further simplified by exponentiating 
the product into a sum, setting the scale as I — 2^f 
and approximating the sum by an integral. We get 

Pr(i? 2 - n < x) cx exp 12^ [ d££ (rjx 2 (11) 
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FIG. 3: Top panel: Pi = Pr(R 2 - fi< x), where R is the 
modulus of the largest eigenvalue, for the open kicked rotator 
at several dimensions, for fi = 0.8. Lower panel: same for 



P 2 = Pr(5f (if 
prediction. 



fj) < x). Solid line is a rescaled RMT 



where we have defined the function 
£(z)=ln' 1 + erf(z) 



(12) 



One interesting question that can be answered at this 
point is, how likely is it that a true gap will exist at /i for 
a finite value of Nl The probability that all eigenvalues 
are smaller than /i is simply given by the exponential in 



for 



(fTTjl with x = 0. It is thus proportional to e 
some constant c. 

Notice that (fTTj) has some similarity with the Tracy- 
Widom distribution [23| of the largest eigenvalue of Gaus- 
sian ensembles of RMT, in the sense that it involves the 
exponential of the integral of a function that satisfies 
a non-linear differential equation, C"(z) — —2z£'(z) — 
(C(z)) 2 . It is not a Painleve transcendent, however. 

Returning to the calculation, let us change variable 
to 5 = i]x 2 and obtain exp |2/i?7 J s C(z)dz\ . Clearly, 
this function does not converge as r\ — > oo. Assuming 
S to be large, we use C{z) w — e -2 /2y/nz and inte- 
grate by parts to get exp |— 2(irie~ s2 /Ay^S 2 ^ . In order 
to have a finite limit, we must have rye - * 5 /5 2 — O(l). 
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This implies that S 2 = y + W(rf), where W is the Lam- 
bert function, which for large rj can be approximated as 
W(rf) ~ log r] — log log J). Therefore, if instead of R we 
consider the variable 
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p = if (R 2 - p) 2 - log x] + log log 77, 



then 



Pr(,o < y) = exp{- 



2^ 



'}• 



(13) 



(14) 



which is a modified Gumbcl function. Therefore, the dis- 
tribution of the slightly awkward variable p (see also [22I ]) 
has a limit as 7/ —¥ 00, but this limit is approached very 
slowly and finite r\ calculations may present significant 
deviations. 

In Figure 3a we see Pr(i? 2 — p < x) for the open kicked 
rotator, at various dimensions for p = 0.8. In Figure 3b 
we introduce the scaled variable rj 2 (R 2 — p) 2 , where rj is 
given by (111). As a result, all curves fall on top of each 
other, indicating that this is the correct scaling. Notice 
that we must not introduce the log rj or log log rj factors 
that appear in (fT3|). as they would spoil the scaling. The 
results agree very well with the function exp{— e~ ax+b }, 
with fitted values of a and b. 

We close with some remarks on resonance eigenfunc- 
tions. These may be depicted in phase space by means 
of their Husimi function, H^(q,p) — \(q, M>)\ , where 



\q,p) is a coherent state. It was shown in [14| that these 
Husimi functions are supported on the backward trapped 
set, the unstable manifold of the repeller. How they are 
distributed on this support depends on their decay rate: 
states with larger T concentrate on the dynamical pre- 
images of the opening, while states with small T concen- 
trate on the repeller. Semiclassically, 



H^(q,p)dqdp a \z n \ 2m (l - \z n \ 2 ), (15) 



where R m is the m-th preimage of the opening. In 
principle, this relation would allow a state whose de- 
cay rate equals the classical decay rate to be uniformly 
distributed, because the area of R m decays like e - ™ 7 . 
Therefore, super-sharp resonances must show an in- 
creased degree of localization above uniformity. Indeed, 
since they can also be seen as super-long-lived states, one 
would expect them to be associated with periodic orbits 
(see for example 24j-|27|). This topic deserves further 
investigations. 

Another line of research to be followed would be to 
investigate a possible relation between the exponents a 
and ft to ray dynamics. In particular, it is not clear 
whether they depend on the Lyapounov exponent and 
whether they are universal or system specific. 
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